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Abstract 

In a recent article the authors obtained a formula which relates explic- 
itly the tail of risk neutral returns with the wing behavior of the Black Sc- 
holes implied volatility smile. In situations where precise tail asymptotics 
are unknown but a moment generating function is available we first estab- 
lish, under easy-to-check conditions, tail asymptoics on logarithmic scale 
as soft applications of standard Tauberian theorems. Such asymptotics 
are enough to make the tail-wing formula work and we so obtain a version 
of Lee's moment formula with the novel guarantee that there is indeed a 
limiting slope when plotting implied variance against log-strike. We ap- 
ply these results to time-changed Levy models and the Heston model. In 
particular, the term-structure of the wings can be analytically understood. 

1 Introduction 



Consider a random variable X whose moment generating function (mgf) M is 
known in closed form, but whose density / (if it exists) and distribution function 
F are, even asymptotically, unknown. For a large class of distributions used for 
modelling (risk-neutral) returns in finance, M is finite only on part of the real 
line. Let us define F = 1 — F and r* as the least upper bound of all real r 
for which M (r) = i?[e rX ] < oo and assume r* £ (0, oo). An easy Chebyshev 
argument gives 

hm sup — = r , (1) 

x— >oc X 

but counter-examples show that the stronger statement 

— \ogF(x) ~ r*x as x — > oo (2) 

may not be true 1 . However, we do expect (2) to be true if the (right) tail of 
the distribution is reasonably behaved. Our interest in such distributions stems 

x We use the standard notation g (%) ~ h (x) = lim g (x) /h (x) = 1 as x —* oo. 
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from the fact that the crude tail asymptotics (2) and the mild integrability 
condition p* = r* — 1 > are enough, via the tail- wing formula [4], to assert 
existence of a limiting slope of Black Scholes implied variance V 2 as function of 
log-strike k. Indeed, in standard notation, reviewed in section 4, one has 

lim V 2 (k)/k = 2-4 ( vV) 2 + P* - P*) • (3) 

k— »OG V / 

Similarly, if q* = sup {q £ R : M (-q) = E[e~ qX ] < oo} £ (0, oo) and the (left) 
tail is reasonably behaved one expects logF(— x) ~ — q*x as x — > oo in which 
case the tail wing formula gives 

lim V 2 (-k) /k = 2 - 4 ( V(q*) 2 +q*- q*) • (4) 

It was already pointed out in [4] that the tail-wing formulae sharpen Lee's 
celebrated moment formulae [9, 8]. In the present context, this amounts to 
having a lim instead of a lim sup 2 . It must be noted that the tail- wing formula 
requires some knowledge of the tails whereas the moment formula is conveniently 
applicable by looking at the mgf (to obtain the critical values r* and — q* ). 

In this paper we develop criteria, checkable by looking a little closer at the 
mgf (near r* and —q*), which will guarantee that (3) rcsp. (4) hold. In view of 
the tail-wing formula the problem is reduced to obtain criteria for (2) resp. its 
left-sided analogue. The proofs rely on Tauberian theorems and, as one expects, 
the monograph [5] is our splendid source. 

The criteria are then fine-tuned to the fashionable class of time-changed 
Levy models [11, 6] and checked explicitly for the examples of Variance Gamma 
under Gamma-OU clock and Normal Inverse Gaussian with CIR clock. We also 
check the criteria for the Heston model. In fact, it appears to us that most (if 
not all) sensible models for stock returns with known mgf and p*,q* £ (0, oo) 
satisfy one of our criteria so that (3) and (4) will hold. 

Finally, we present some numerical results. The asymptotic regime becomes 
visible for remarkably low log-strikes which underlines the practical value of 
moment - and tail-wing formulae. 

2 Background in Regular Variation 

2.1 Asymptotic inversion 

If / = / (a;) is defined and locally bounded on [X, oo), and tends to oo as x — > oo 
then the generalized inverse 

(x) := inf {y £ [X, oo) : / (y) > x} 

2 Remark that, at least when p* > 0, the moment formula is in fact recovered from the 
tail-wing formula and (1). 
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is defined on [/(A),oo) and is monotone increasing to oo. This applies in 
particular to / e R a with a > and Thm 1.5.12 in [5] asserts that /*~ e R\/ a 
and 

/ (/^~ {x)) ~ (/ (a;)) - z as a; -> oo. 

Given / one can often compute f*~ (up the asymptotic equivalence) in terms 
of the Bruijn conjugate of slowly varying functions (Prop. 1.5.15, Section 5.2. 
and Appendix 5 in [5]). 

2.2 Smooth Variation 

A positive function g defined in some neighbourhood of oo varies smoothly with 
index a, g € SR a , iff h (x) := log (g (e x )) is C°° and 

h! (x) — > a, M™) (x) — ^ for n = 2, 3, ... asu oo. 

Theorem 1 (Smooth Variation Theorem, Thm 1.8.2 in [5]) If f e R a 

then there exist /, g S'i? ct; i = 1,2, loitA /i ~ / 2 and /i < / < / 2 on some 
neighbourhood of oo. 

When a > we can assume that /i and / 2 are strictly increasing in some 
neighbourhood of oo. In fact, we have 

Proposition 2 Let a > and g g S'i? ct . T/ien g is strictly increasing in some 
neighbourhood of oo and g' g SR a -\. 

Proof. By definition of S'-Rq, 

9 1 

— log (g (e x )) = — r— V (e x ) e I ^a>0asi^oo. 
ox g(e x ) 

This shows that, in some neighbourhood of oo, g' is strictly positive which 
implies that g is strictly increasing. From Prop 1.8.1 in [5], g' = \g'\ g SR a -\. 
■ 

Remark 3 In the situation of the last Proposition we have lim^oo g (x) = oo 
and hence, in some neighbourhood of oo , g has a genuine inverse g^ 1 which 
coincides with the generalized inverse g*~ . 

2.3 Exponential Tauberian Theory 

Theorem 4 (Kohlbecker's Theorem, Thm 4.12.1 and Cor 4.12.6 in [5]) 

Let U be a non- decreasing right- continuous function on R with U (x) = for all 
x < 0. Set 

N(X) := [ e- x / x dU{x), A > 0. 

J[0,oo) 

Let a > 1 and \ S R a /{a-i)- Then 

\ogN(X) ~ (a - 1) x (A) /A as A -> oo 
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iff 

log fi [0, x] ~ axjx~ (x) as x — > oo. 

Theorem 5 (Karamata's Tauberian Theorem, Thm 1.7.1 in [5]) LetU 
be a non- decreasing right- continuous function on R with U (x) = /or all x < 0. 
IflGRo an d c > 0, p > 0, £/ie following are equivalent: 

U (x) ~ cx p Z (a;) /r (1 + p) asi^oo 

t/(s)= / e- s:E dJ7(a;) - ca-"l(l/a) as s -» + . 
Jo 

(When c = the asymptotic relations are interpreted in the sense that U (x) = 
o (x p l (x j) and similar for U.) 

Theorem 6 (Bingham's Lemma, Thm 4.12.10 in [5]) Let f e R a with 
a > such that that e~f is locally integrable at +oo. Then 

/•OO 

-log / e - f Mdy ~ f(x) . 

J X 

3 Moment generating functions and log-tails 

Let F be a finite Borel measure on R, identified with its (bounded, non- 
decreasing, right-continuous) distributions function, F (x) = F ((— oo, a;]). Its 
mgf is defined as 

M (s) := J e sx dF (x) . 
We define the critical exponents q* and r* via 

— q* = inf {s : M (s) < oo} , r* = sup {s : M (s) < oo} 
and make the standing assumption that 

r*,q* E (0,oo). 

In this section, we develop criteria which will imply the asymptotic relations 

\ogF ((— oo, —a;]) ~ —q*x, \ogF {fx, oo)) ~ —r*x asu oo. 

The assumption in the following Criterion I is simply that some derivative 
of the mgf (at the critical exponent ) blows up in a regularly varying way. 

Theorem 7 (Criterion I) Let F be a bounded non- decreasing right- continuous 
function on R and define M = M (s) , q* and r* as above. 

(i) If for some n> 0, M (n) (—q* + s) — s^ p li(l/s) for some p > 0, h £ R as 
s — > 0+ then 

logF ((— oo, —x\) ~ — g*a; 

fiij ///or some n > 0, M(") (r* - s) ~ s~ p li(l/s) for some p > 0, li E R as 
s — > 0+ i/ierc 

log F ((#, oo)) ~ — r*x. 
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Proof. Let us focus on case (ii), noting that case (i) is similar. We first discuss 
n = 0. The idea is an Escher-type change of measure followed by an application 
of Karamata's Tauberian Theorem. We define a new measure U on [0, oo) by a 
change-of-measure designed to get rid of the exponential decay, 

dU (x) := exp (r*x) dF (x) . 

We identify U with its non-decreasing right-continuous distribution function 
x i— > U ([0, x}). The Laplace transform of U is given by 

/•oo />oo r-0 

U{s)= e~ sx dU(x) = e^'- s)x dF{x)=M(r* - s)- e^*-'>dF (x) 

JO JO J -oo 

so that 

r° 

U(s)-M(r*-s) < / e ( - r "- s '> x dF(x)<F(0)-F(-oo)<2\\F\\ oct <oo. 

J —oo 

Since M (r* — s) goes to oo as s — > 0+ and we see that U (s) ~ M (r* — s) so 
that U e R p as s -> 0. Hence, there exists I £ R so that {/ (s) = (l/s) p / (1/s) 
and Karamata's Tauberian theorem tells us that U e i? p , namely 

C/ (a:) - a; p / (a;) /T (1 + p) = x p l' (x) as x -» oo 

where I' £ Rq. Going back to the right-tail of F, we have for a; > 



F((a;,oo)) 



f dF(y)= f exp(-r*y)dU(y). 



We first assume that U € SR p . Under this assumption U is smooth with 
derivative u = U' e SR p -i and we can write 

u (y) = y"- 1 /" (y) with /" e R - 

Then 



F((ar,oo))= / exp (-r*y) y^H" (y) dy 

J (x,oo) 

= / exp [-r*y + (p - 1) log y + log I" (y)] dy. 

J (x,oo) 



' (a:,oo) 

Since — [—r*y + (p — 1) logy + log I" (y)] ~ r*y e i?i as y — > oo we can use 
Bingham's lemma to obtain 

-logF((a;,oo)) = -log / exp [-r*y + (p - 1) logy + log I" (y)] dU (y) - r*y. 

J (x,Cso) 

(5) 
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We now deal with the general case of non-decreasing U € R p . From the Smooth 
Variation Theorem and Proposition 2 we can find U-,U + G SR P , strictly in- 
creasing in a neighbourhood of oo, so that 

[/_<[/< U+ and U- ~ U ~ f7+. 

Below we use the change of variable z = U (y) and w = l^ 1 (z). Noting that 
U+ 1 < < UZ 1 and using changc-of- variable formulae, as found in [10, p7-9] 
for instance, we have 

F((x,oo)) = f exp{-r*y)dU(y) 

J (x.oo) 

= [ exp(-r*[/^ {z))dz 



< 



[ exp(-r*C/+ 1 (z))rfz 
= / exp (— r*w) dU + (w) . 

J(U-\U(x)), X ) 

Similar to the derivation of (5), Bingham's lemma leads to 

- log f exp (-r*w) dU+ (w) - r*U7 1 (U (x)) . 

J(U^(U(x)),^) 

Noting that UI 1 is non-decreasing, U+ 1 (U (x)) < U+ 1 (U + (x)) = x so that 3 

— logF ([x, oo)) < r*x 

The same argument gives the lower bound — logF ((a;, oo)) > r*x and we con- 
clude that 

— logF ((a;, oo)) ~ r*x. 

We now show how n > follows from n = 0. Define V on [0, oo) by 

dV (x) := x n dF (x) . 

Clearly, V induces a non-decreasing, right continuous distribution on R, V (x) := 
V ([0, x]) for x > and V (x) = for x < 0. The distribution function V (x) is 
also bounded since 

/■OO 

/ (x) < oo 

Jo 

which follows a forteriori from the standing assumption of exponential moments. 
We will write V(x) for V(x, oo). 

3 By g < h we mean limsup / (x) jg (x) < 1 as x — + oo. 



G 



Note that V has a mgf My (s), finite at least for s€ (0, r*), given by 
My (s) = f e sx dV (x) = f x n e sx dF 




o 



x n e sx dF + C = tfl") (s) + C 



where' 



,4 




J —oo J —oo 



By assumption, M ^ is regularly varying with index p at r* and it follows that, 
as s — > 0+, 



We now use the "n = 0" result on the distribution function V resp. its mgf My 
and obtain 

— log V ([x, oo)) = — log V (x) ~ r*x e R\ 
Assume first that — log V (x) G SRi . Then V has a density V = v and 

v(x) = d x {V(oo) - V{x)) = -V{x)d x (log? (a;)) - r*V (x) as x -> oo 

since functions in Si?! are stable under differentiation in the sense that d x (— log F (x)) ~ 
9a; (r*x) = r*. In particular, we have \ogv (x) ~ log V" (x) ~ — r*x. After these 
preparations we can write 



and Bingham's lemma implies that logF((x,oo)) ~ — r*x. The general case 
of log V (x) G Ri follows by a smooth variation and comparison argument as 
earlier. ■ 

The next criterion deals with exponential blow-up of M at its critical values. 

Theorem 8 (Criterion II) Let F, M ', q* , r* be as above. 

(i) If log M (—q* + s) — s~ p li(l/s) for some p > 0, h e Rq as s — » 0+ then 



(ii) If log M (r* — s) ~ s p l\{l/ s) for some p > 0, h G Rq as s — ► 0+ then 



M v {r* -s) = M(") (r* - a) + O (1) - a- p h(l/a). 





exp[logv(y) -nlogy]dy 



logF((-oo,-x]) 



logF ((x, oo)) ~ — r*x. 



4 One could do without the assumption J, ^ \x\dF (which follows a fortcriori from the 
standing assumption q* > 0). Finitcncss of F on (— oo,0) is enough. 
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Proof. As for Criterion I, the idea is an Escher-type change of measure followed 
by a suitable Tauberian theorem; in the present case we need Kohlbecker's 
Theorem. Let us focus on case (ii), noting that case (i) is similar. A new 
measure U on [0, oo) is defined by 

dU (x) := exp (r*x) dF (x) . 

We identify U with its non-decreasing right-continuous distribution function 
x \— ► U ([0, a;]) and define the transform 

/■OO P GO [>0 

N (A) = / e- x ' x dU (x) = / e^'-W'dF (x) = M (r* — 1/A)- / e^^'^dF (x) 

JO JO J -oo 

so that 

,o 

\N (A) - M(r* - 1/A)| < / -W'dF (x) < F (0)-F (-oo) < 2 < oo. 

J — OO 

Thus, 

N (A) = M(r*- 1/A) + 0(1) as A ^ oo 

and, in particular, since liniA^oo logM (r* — 1/A) = limA^oo M (r* — 1/A) = oo 
from the assumption (ii) we see that 

log A (A) - log M(r* - 1/A) - A p /i(A) £ R P as A -> oo. 

Define a e (1, oo) as the unique solution to p + 1 = a/ (a — 1) and note 

X(A) := — ^ - log TV (A) £ iJp+i = J? Q /( a _i). 

Using that € R( a -i)/a = R\-i/ a i Kohlbecker's Tauberian Theorem tells us 
that 

\ogU ([0, x]) = log/7 (x) ~ ax/x^ (x) G i?i/ Q as x — > oo. 

In particular, there exists I € i?o so that log [/ (x) = ax 1 /"? (x). We first assume 
that log/7 £ SRi/ a . Then /7 has a density u (.) £ 5i? 1 / ct _ 1 and 

u (x) = U (x) 9 X (log U(x))~U (x) x 1 /"- 1 / (x) . 

In particular, 

logu (x) ~ logf/ (x) e -Ri/ Q as x — > oo. 
Now, y i ► r*y G i?i dominates i?i/ Q (since 1/a < 1) in the sense that 

r*y- logu(y) <~ r*y. 

Thus, from 

F((x,oo))= dF(y) = exp(-r*y)u(y)dy 

J (x.oc) J [x,oo) 

exp [-r*y + \ogu(y)] 



L 



(x,oo) 
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and Bingham's lemma we deduce that 

— logF ((x, oo)) ~ r*x. 

The general case, log U e Ri/ a , is handled via smooth variation as earlier. 
Namely, we can find smooth minorizing and majorizing functions for log U, say 
G_ and G+ . After defining U± = exp G± we have 

log U- ~ log U ~ log C/+ and U-<U <U+. 

Then, exactly as in the last step of the proof of Criterion I, 

F((x,oo))= / exp (-r*y) dU (y) < / exp (-r*w) dl/ + (w) 

J(x,oo) J(U^ 1 (U(x)),oo) 

and from Bingham's lemma, 

- \ogF {(x, oo)) < r*£/+ 1 (U (x)) ~ r*x. 
Similarly, — \ogF ((x, oo)) > r*x and the proof is finished. ■ 

4 Application to Smile Asymptotics 

We start with a few recalls to settle the notation. The normalized price of a 
Black-Scholes call with log-strike k is given by 

c BS (k,a) = $(d 1 )-e k $(d 2 ) 

with c?i j2 (fc) = —k/a± a/2. If one models risk-neutral returns with a distribu- 
tion function F, the implied volatility is the (unique) value V (k) so that 

/>OC 

c BS (k, V (k)) = (e x - e k ) dF (x) =: c (k) . 

Jk 

Set ip [x] = 2 — 4 [V x 2 + x — x) and recall F = 1 — F. The following is a special 
case of the tail- wing formula [4]. 

Theorem 9 Assume that — logF (— k) jk ~ q* for some q* € (0, oo) . Then 
V(-kf/k - ^ [- \ogF(-k) /k] - V (q*) ■ 

Similarly, assume that — logF (k) jk ~ p* + 1 /or some G (0, oo). Tften 

F(fc) 2 /fc-^[-l-logF(fc)/fc] ~ 

As earlier, let M (s) = J exp (sx) di 7, (x) denote the mgf of risk-neutral re- 
turns and now define the critical exponents r*and — q* exactly as in the begin- 
ning of the last section 3. Combining the results therein with Theorem 9 we 
obtain 
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Theorem 10 If q* € (0,oo) and M satisfies part (i) of Criteria I or II then 

V(-k) 2 /k~iP(q*) ask^oo. 

Similarly, if r* = p* + 1 G (l,oo) and M satisfies part (ii) of Criteria I or II 
then 

V(k) 2 /k~ip{p*) ask^oo. 

5 First Examples 

The examples discussed in this section model risk-neutral log-price by Levy 
processes and there is no loss of generality to focus on unit time. 5 

5.1 Criterion I with n=0: the Variance Gamma Model 



The Variance Gamma model VG = VG (m, g, C) has mgf 
M(s) 



c 

gm \ I gm 



gm + (m — g)s — s 2 J \ (m — s)(s + g) 

The critical exponents are obviously given by r* = m and q* = g. Focusing on 
the first, we have 

M(r* - s) ~ ( ) s- c as s -> 0+ 

which shows the Criterion I is satisfied with n = 0. Theorem 10 now identifies 
the asymptotic slope of the implied variance to be if> (r* — 1) = ip(m— 1). 
Similarly, the left slope is seen to be tp (q*) — ip (g). We remark that [1] contains 
tail estimates for VG which lead, via the tail-wing formula, to the same result. 

5.2 Criterion I with n>0: the Normal Inverse Gaussian 
Model 

The Normal Inverse Gaussian Model NIG = NIG (a, [3, /i, 5) has mgf given by 

M (s) = exp | S jVa 2 - (3 2 - \J a 2 - ((3 + s) 2 | + ^sj . 

By looking at the cndpoints of the strip of analyticity the critical exponents are 
immediately seen to be r* = a — (3, q* = a + (3 and we focus again on the first. 
While M(s) converges to the finite constant M(r*) as s — > r*— we have 

M'(s)/M (s) = (25((3 + s)[a 2 -([3 + s) 2 ]- 1 ' 2 + M ) 

and M'(r* - s) - 2<W2aV 1/2 : M(r*) as s + . 



5 In fact, Levy models that satisfy one of our criteria have no term structure of implied 
variance slopes. 
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We see that Criterion I is satisfied with n = 1 and Theorem 10 gives the 
asymptotic slope tp (r* — 1) = ip (a — (i — 1). Similarly, the left slope is seen to 
be V (q*) — 4> (& + P)- We remark that the same slopes were computed in [4] 
using the tail-wing formula and explicitly known density asymptotics for NIG. 

5.3 Criterion II: the Double Exponential Model 

The double exponential model DE — DE (a, /x, X,p, q, r/i, 772) has mgf 



and we see that Criterion II is satisfied. As above, this implies asympotic slopes 
tp (r* — 1) = ip (r/i — 1) on the right and ip (772) on the left. 

6 Time changed Levy processes 

We now discuss how to apply our results to time changed Levy processes [11, 
12, 6]. To do this, we only need to check that the moment generating function 
of the marginals of the process, satisfies one of the three criteria. 

To this end, consider a Levy process L = L(t) described through its cumu- 
lant generating function (cgf) Kl at time 1, that is, 



and an independent random clock T = T (oj) > with cgf Kt- It follows that 
the mgf of L o T is given by 



Therefore, in order to apply our Theorem 10 to time-changed Levy models 
we need to check if M = exp [Kt(Kl(-)} satisfies criterion I or II so that 
— \ogF(x)/x tends to a positive constant. Here, as earlier, F denotes the 
distribution function of M and F = 1 — F. The following theorem gives suffi- 
cient conditions for this in terms of Kt and Kl. We shall write Mt = exp (Kt) 
and Ml = exp (Kl). For brevity, we only discuss the right tail 6 and set 



°In fact, the elegant change-of-measure argument in Lee [9] allows a formal reduction of 
the left tail behaviour to the right tail behaviour. 




Clearly, r* — rji and as s — > 0+ 




K L (v) = log E [exp (vL x )] 



M(v) = E [E (e vLT \T)] = E e KL ^ T = exp [K T (K L (v)] . 



PL = sup {s : M L (s) < 00} , p T = sup {s : M T (s) < 00} . 
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Theorem 11 With notation as above, assuming pl^Pt > 0, we have: 

(1.1) If Kl(jj) = pt for some p <G [0,Pl) and Mt satisfies either criterion then 

logF(x) ~ — px. 

(1.2) If Kl(j>) = Pt for p = pl and Mt,Ml satisfy either criterion then 

logF(x) ~ — px. 

(ii) If Kl{p) < Pt for all p e [0,pl] and Ml satisfies either criterion then 

logF(x) ~ -Plx. 

Remark 12 It is worth noting that there cannot be more than one solution 
to Kl(p) = pt- To see this, take any v such that v > and Kl (v) > (any 
solution to Kl(p) — pt > will satisfy this!) From Ml = exp Kl it follows that 
Ml (0) = 1 and Ml (v) > 1. By convexity of Ml (•) it is easy to see that M' L (v) 
is strictly positive and the same is true for K' L (v) = M' L (v) /Ml (v). It follows 
that w > v => K L (w) > K L (v) > and the set of all v > : K L (v) > 
is connected and Kl restricted to this set is strictly increasing. This shows that 
there is at most one solution to Kl{p) — Pt- 

Proof, (i.l) Noting that p > let as first assume that Mt satisfies criterion I 
(at Kl{p) — Pt with some n > ) so that for some p > and I £ Rq, 

M T n ^ (u) <~ (pt — u)~ p I {(pt — u) 1 ^ as u t pt- 

From M = M T ° K L we have M' = M T {K L )K' L and, by iteration, M<"' 
equals M [ T n) (K L ) {K' L ) n plus a polynomial in M T (•) , ...,M { T n ~ 1] (•) which re- 
mains bounded when the argument approaches pt- Noting that K' L (p) > (see 
remark above) we absorb the factor [K' L (p)] n into the slowly varying function 
and see that 

M (n) (w) - (p T - K L (v))~ p l((p T - KL^y 1 ) for p as above and some I e Ro 

as Kl (v) tends to pt which follows from v f p. Using analyticity of Kl in 
(0,pl) and K' L (p) ^ it is clear that pt — Kl(v) ~ K' L (p)(p — v) as v | P and 
so 

M {n \p - v) ~ K' L (p)- p v- p l(l/v) as v -» + . 

This shows that M satisfies Criterion I (with the same n as Mt). A similar 
argument shows that M satisfies Criterion II if Mt does. Either way, the assert 
tail behaviour of logF follows. 

(i.2) The (unlikely!) case Kl(pl) = Pt involves similar ideas and is left to the 
reader. 

(ii) We now assume that 

sup K L {p) < Pt < oo. 
pe[o,Pi] 
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and Ml satisfies either criterion (at pl)- Since Ml = cx^Kl stays bounded as 
its argument approaches the critical value pl it is clear that Ml cannot satisfy 
criterion II or criterion I with n = and there must exist a smallest integer n 
such that 

M L n \pL — x) ~ x~ p l(x) as x — ► pl 
for some p > and I E Rq. We note that 

M™(v) = (^( V )4(ir L («)) + /(«)) ex V {K T {K L {v))) 

where f(v) is a polynomial function of the first (n — 1) derivatives of Kl and 
the first n derivatives of Kt evaluated at Kl(v), which are all bounded for 
< v < pl- Noting that positivity of T implies M' T > and hence K' T > we 
see that as v | Pl 

M^ n \v) ~ K^^K^KLiPL^M ( PL ) . 

Applying this to Kt (x) = x leads immediately to 

K L n \v) ~ M L n \v)/M L (v) ~ i^(i)/M L (p L ). 

as w | and so M satisfies criterion I. 
■ 

We now discuss examples to which the above analysis is applicable. For all 
examples we plot the total variance smile 7 for several maturities and compare 
with straight lines 8 of correct slope as predicted by Theorem 10. All plots are 
based on the calibrations obtained in [12]. This is also where the reader can 
find more details about the respective model parameters. 



6.1 Variance Gamma with Gamma-OU time change 

We will consider the Variance Gamma process with a Gamma-Ornstein-Uhlenbeck 
time change and refer to [12] for details. From earlier, the Variance Gamma pro- 
cess has cumulant generating function 

K L {v) = Clog (- 9 ™ ) for v e {-g, m) 

\(m - v)(v + g) J 

We note that Kl ([0, m]) = [0, oo] so that pl = oo. The Gamma-Ornstein- 
Uhlenbeck clock T = T{uj, t) has cgf 

K T (v) = vy \- 1 (l-e- xt )+ A " 



v-Xb 



bl0 ^ b- V X-^l-e-^) l- Vt 



We need to examine how this function behaves around the endpoint of its strip 
of regularity. At first glance, it appears that the function tends to infinity as 

7 That is, V 2 (k, t) = a 2 (k, t) t. 

8 These lines have been parallel-shifted so that they are easier to compare with the actual 
smile. 
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v t A6, because of the ^jz^j; term. However, upon closer examination, we can 
see that this is in fact a removable singularity, and the term of interest to us 
is the log(...) term. This term tends to infinity as v — > A6(l — e~ At ) _1 =: pr- 
After some simple algebra, we see that 

PT A PT ~ xb f_ „_ A n PT^Aa 



exp<^p T j/oA (1-e ) — ^ as v f Pt 



Pt — v J ( pr — Afo 

Therefore, exp^y) satisfies Criterion I with n = and part (i.l) of Theorem 
11 shows that M does too and so that logF(x) ~ —px where p is determined 
by the equation 

K L (p)=p T = M(l-e- Xt )- 1 
and can be calculated explicitly, 

_ to - g + y/(m - g) 2 + 4gm(l - exp(-A6/C(l - e~ A *)) 




-2-1 12 

Figure 1: VG with Gamma-OU time change. Parameters from [12]. Total 
implied variance and slopes for three maturities t = 0.4, 0.9 and 1.3 years. 



6.2 Normal Inverse Gaussian with CIR time change 

The cgf of the Cox-Ingersoll-Ross (CIR) clock T — T (u>, t) is given by 

K T (v) = K 2 77VA 2 + 2y w/(K + 7Coth(7i/2))- ^log[sinh7i/2(coth7t/2 + -) 

A 7 
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where 

7 = \/ n 2 — 2X 2 v. 

This clearly tends to infinity as I(v) = n + j(v) coth(7(w)t/2) — » 0, and we can 
define px as solution to the equation I (px) = 0. Using l'Hopital's rule, it is 
easy to check that 

Pt-v 

K + 7(w)coth(7(w)t/2) 

tends to a constant as v —* , and so 2yow/ (k+7 coth(7*/2) is regularly varying 
of index 1 as a function of (pt — v) 1 . ft is clear that this is the dominant term 
in this limit, and so Mr = exp (Kt) satisfies criterion II (at pr)- From earlier, 
the NIG cgf is 9 

K L (v) = -S(^a 2 - (fi + v) 2 - ^a 2 - f3 2 ) for v < a - (3 

from which we see that pl = a — (3 > and 

sup Kl(v) = 5\/a 2 — (3 2 . 

ve[0,a-/3] 

Therefore, the behavior of M on the edge of the strip of analyticity, and the 
location of the critical value, will depend on whether this supremum is more or 
less than pr; if it is less than pt, the latter is never reached. Recalling that 
exp (Kl) satisfies Criterion I with n = 1, we apply part (ii) of Theorem 11 and 
obtain 

— log^(x) <~ plx = (a — [i)x. 

Otherwise, there exists p e (0, a — (5} such that Ki,(jp) = pt, for some p < a — {3, 
and since Mt was seen to satisfy one of the criteria (to be precise: Criterion II) 
we can apply part (i) of Theorem 11 and obtain 

— logF(x) ~ px. 

In particular, we see that for all possible parameters in the NIG-CIR model the 
formula (2) holds true. Smile-asymptotics are now an immediate consequence 
from Theorem 9. 

6.3 The Heston Model 

The Heston model is a stochastic volatility model defined by the following 
stochastic differential equations: 

^ = ^IdW} 
dv t = k(t] — v t )dt + v t dW t 

9 Following [12] we take fi = here. 
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Figure 2: NIG with CIR time change. Parameters from [12]. Total implied 
variance and slopes for three maturities t — 0.4, 0.9 and 1.3 years. Observe that 
the the lines with correct slope do not perfectly line up with the smile which is 
not a contradiction to the result that that V 2 (k) jk converges to a constant. 



where d(W^, Wf\ — pdt is the correlation of the two Brownian motions. log St 
therefore has the distribution of a Brownian motion with drift —1/2 evaluated 
at a random time T {uj, t) = J Q v s ds with the distribution of an integrated CIR 
process, as in the previous example. When p = 0, the Levy process L = W 1 
and T are independent and we can apply the same analysis as above. Namely, 
the cgf of the Brownian motion with drift speed —1/2 at time 1 is 

K L {v) = {v 2 -v)/2, 

so that pl — oo, and Mt = exp (Kt) satisfies Criterion II hence, by part (i) of 
Theorem 11, 

logF(x) ~ —px 

where p is determined by the equation Kl (p) = pj>. When p < 0, we can 
analyze the mgf of log St directly, and we can apply the same reasoning as 
for the mgf of the CIR process, to deduce that criterion II is satisfied. The 
distribution function for the Heston returns hence satisfies \ogF(x) ~ — px 
where p is solution to, see [3], 



(« - pvff) + {0 2 (v 2 -«)-(«_ pvO) 2 ) 1 ' 2 cot{{e 2 {v 2 -v)-(k- pvQfflHjl} 



v—p 



When p > 0, which is of little practical importance (at least in equity mar- 
kets), the mgf may explode at a different point, see [3], but criterion II will still 
be satisfied. 
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Figure 3: Heston Model. Parameters from [12]. Total implied variance and 
slopes for three maturities t = 0.4, 0.9 and 1.3 years. 
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